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Abstract 

It is shown that the geometric measure of entanglement of a pure 
multipartite state satisfies a polynomial equation, generalising the 
characteristic equation of the matrix of coefficients of a bipartite state. 
The equation is solved for a class of three-qubit states. 

1 Introduction 

A pure state of a composite quantum system is entangled if it is not a product 
of states of the individual parts of the system. Thus one way of measuring the 
entanglement of the state (though there are many other independent types of 
entanglement [9]) is to measure how different the state is from any product 
state. If we consider only pure states, this can be expressed in geometrical 
terms as the distance from the closest pure product state. The space of pure 
states of a system described by a Hilbert space Ti is the projective space 

of one dimensional subspaces of 7i, and the closeness of two pure states 
-01 and -02 is measured by the modulus of the inner product |(V^i|V'2)| where 

and \ip2) are normalised vectors representing the states; this gives rise 
to a standard metric, the Fubini-Study metric [13], on the projective space 
P7Y. In the multipartite situation the Hilbert space is a tensor product 
Ti = Til ® ■ ■ ■ ® Tin, and the geometric measure of entanglement [12, 11] or 
Groverian measure [1] of a pure state E H is a. decreasing function of 



its Fubini-Study distance g from the set of product states |0i) • • • with 

5(|*)) = max|(*|(|0i)---|0„))|2 

where the maximum is taken over all single-system normalised state vectors 
|0i) G Til, • • • 5 G Ti„, and it is understood that |\E') is normalised. 

The function g is known as the injective tensor norm on the tensor prod- 
uct space Til ® ■ ■ ■ ® Ti-n- In the bipartite case (n = 2), ^id^^)) is the largest 
coefficient in the Schmidt decomposition of |^) (equivalently, the largest 
eigenvalue of the density operator |^)(^|, or the largest singular value of the 
matrix of coefficients of |^) with respect to a product basis of Tii (S> ^^2); 
for general n it is the first coefficient in the multipartite generalisation of 
the Schmidt decomposition proposed in [2]. It is determined by equations 
which constitute a nonlinear generalisation of the equations for the singu- 
lar values and singular vectors of a matrix, the matrix being replaced by a 
"hypermatrix" (an n-dimensional array of numbers, here the coefficients of 
the given state with respect to a product basis). For these reasons, gd"^)) 
is also known as the "maximum eigenvalue" (more correctly, the "maximum 
singular value") of the state 1^^). 

The definition of gi as a maximum makes it susceptible to numerical cal- 
culation, but it is not so easy to treat it analytically. It has been determined 
in certain special cases [12, 11, 14, 10, 5, 8] but no general method of calcu- 
lation has been proposed. In this paper we will give such a method, showing 
how the problem can be reduced in principle to the solution of a polynomial 
equation. We apply the method to three-qubit states. However, the polyno- 
mial equation in question has a very high degree, and even for three qubits 
the general case remains to some extent intractable. For a special class of 
three-qubit states, we obtain a complete solution, reproducing the results 
obtained by other authors [10] by methods which are specific to qubits. 

The equations for the singular values and singular values of a hypermatrix 
are of mathematical interest in their own right. There does not seem to be 
a general theory available. We conclude by comparing the results for the 
particular case studied here with the familiar theorems for the matrix case. 
It appears that few of these theorems will have straightforward analogues in 
the general case. 
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2 The nonlinear singular- value equations 

Given |\E') G ?ii ■ ■ ■ 7i„, we want to find \4>k) G Ti-k, k = 1, . . . ,n, which 
maximise the function 

subject to {(j)k\4>k) — 1- Let le^*^^) {i = 1, . . . ,dr — diuiHr) be an orthonormal 
basis ior Tir {r — 1, . . . ,n), and write 

The problem becomes: 

I 2 

Maximise /(u) = laj^.^i^-wj^^^ . . .ufj 
subject to ^ \'^tr\^ ~^ (r = 1, . . . n). 

ir 

where summation is understood over repeated lower indices ii, . . . , i„. Intro- 
ducing Lagrange multipliers A,., we are led to the equations [2, 12] 



where iif^ denotes that n-'-* is to be omitted, and the overline denotes the 
complex conjugate. Multiplying by u^[^ and using the constraints, wc find 
that the have a common value A, which is the required extreme value of 
/(u). We will say that A is a singular value of the hypermatrix ai^_„i^, for 
the following reason. 

In the simple case n = 2 the coefficients a;^;^ form a di x ^2 matrix A, the 
unknowns m-J-* and uf^ form a (ii-component vector u and a ci2-component 
vector V, and the equations can be written 

Av = An, (2.2) 
u^^ = Avt (2.3) 
or A^vi = Av (2.4) 

where the dagger denotes the hermitian conjugate. Multiplying the first 
equation by A^ and using the second shows that v is an eigenvalue of A^A 
with eigenvalue A^, i.e. A is a singular value of A. It satisfies the equation 

det(^U - A'7) = 0. 
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Our aim is to find the corresponding equation in the general case. 

If A = there is an established general theory for the equations (2.1). 
They can be written in terms of the multilinear function a of n vector vari- 
ables defined by 

as 

a = 0, V^a = (r = l,...,n) (2.5) 

where denotes the gradient with respect to the vector variable u^. A 
solution u^^\ . . . ,u(") of these equations is a critical point of the function 
a. The number of independent equations is one more than the number of 
variables, so the variables can be eliminated to give a polynomial equation 
in the coefficients aji...i„: 

Theorem (Cayley)[4] There is a function hdet(a), polynomial in the 
coefficients aj^...i„, such that the equations (2.5) have a solution with all Uj 
non-zero if and only if hdet(Q;) = 0. 

The function hdet(a) is called the hyperdeterminant of the hyperma- 
trix aij^,,,i„- It is a special case of the discriminant [4] of a function of sev- 
eral variables: the discriminant of a homogeneous polynomial function 
f{xi, . . . , xn) is a polynomial in the coefficients of / with the property that 
Aj = if and only if there is a point x at which / vanishes together with all 
its derivatives. The hyperdeterminant has appeared elsewhere in the theory 
of multipartite entanglement [6, 7]; for n = 3 and di = d2 = = 2, its 
squared modulus is equal to the 3-tangle [3] of the 3-qubit state defined by 
the coefficients aijk. For n = 2 and di — d2, the hyperdeterminant reduces to 
the familiar determinant of the square matrix Oi^i^. In this case, as we have 
already noted, the determinant function also gives the condition for the equa- 
tions (2.1) to have a solution for any A (not just A = 0); the characteristic 
equation of a square matrix, and the singular- value equation of a rectangular 
matrix, are both given by determinants. In general, however, the character- 
istic equation of a tensor is not so simply related to the hyperdeterminant 
function, and the theory needs to be extended. 

3 Generalising the Characteristic Polynomial 

In this section we shall derive a polynomial expression in the coefficients of 
a multipartite pure state, the roots of which include the geometric measure 
of entanglement of the state. This is directly analogous to the characteristic 
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polynomial of a square matrix, which can be used to determine the Schmidt 
coefficients of an ordinary bipartite vector. 

Our method is to consider the function a(u*^^\ . . . , u*^")) as a function 
of two of the variables with a matrix of coefficients which depend on the 
other variables. First we shall prove a lemma concerning such matrix- valued 
functions of several variables. 

Theorem 3.1. Let A{z) be a p x q matrix function of N complex variables 
z = {zi, . . . , zn) which depends only on the variables z^ and not on their 
complex conjugates Zk, i-c the entries of A{z) are holomorphic functions of 
the Zk'; let f3(z) be a (non-holomorphic) real-valued function of z £ C''^; and 
let X be a positive real number. Then the real-valued function 

det[A(z)U(z) - A'/3(z)7] 

vanishes along with its first derivatives with respect to the real and imaginary 
parts of Zk, at a point z with /3(z) = 1, if and only if \ is a singular value of 
A{z) and 

xT^y = (3 1) 

dzk dzk 

where x G and y G C"^ are respectively left and right singular vectors of 
A{z) corresponding to the singular value A; 

A{z)y = Ax, (3.2) 
xT^(z) = Ayt, (3.3) 
xtx = yV = 1. 



Note: in eq. (3.1) /3 is to be treated as a function of independent variables 

, _ d(3 I f dp .dp\ , . , 

z and z (i.e. ^ = ^ ^tt- where Zk^Xk + lyk)- 

ozk I \oXk Oyk J 

Proof. Let the singular values of A{z) be (7i(z), . . . , cr^(z). Then we can write 
the determinant as 

P(z, A) := det[A(z)U(z) - \^(3{z)\ = H " ^''^(^)) " 

i 

Calculating the derivative we have 

aP(z,A) sr^(d[af\ ^,dp 



dzk 



assuming /3(z) = 1. Now P(z, A) = is equivalent to (7i(z)^ = A^, for some i, 
so (Jj(z) = A. Suppose i = 1. If cri(z)^ is a degenerate eigenvalue of A^zY A{z) 
then the derivative will vanish, as the product term on the right hand side 
will always contain a zero term. Thus we shall assume that the eigenvalue is 
non-degenerate. In that case we have 

dP{z,X) _ fd[al] ^,df3 



dzk 



The product term is non-zero, so 

dzk dzk dzk 

Now consider x and y as functions of z and z defined by the singular- value 
equations 

A{z)y = (7i(z)x, 
xTA(z) = <7i(z)yt, 
xtx = yV = 1- 



These give 
and therefore 



M = yt.4.|ly + S^AKAj + yt^,^|y (3.5) 

dzfe a^fe dzk ozk 

^criyi'—y + ai[—y + y^—] (3.6) 

since y'l'y = 1 is constant. The theorem now follows from (3.4) (since P is 
real, dP/dzk = implies the vanishing of derivatives with respect to the real 
and imaginary parts of Zk). □ 

Now take p = di, q = d2 and N = ds + h with x = u*^"^'', y = u^'^\ 

z = (u(3),...,u('^)), and 



/3(z) = ||u(3)||^..||u(")|p 
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where ||u|p = u^u. Then the equations (3.1)-(3.3) become the nonhnear 
singular- value equations (2.1). Thus the nonlinear singular- value problem for 
n vector variables has been reduced to the problem of finding critical points of 
a function of n — 2 vector variables. The function dct[A(z)M(z) — A^/3(z)/] 
is a homogeneous polynomial function of the real and imaginary parts of 
the n — 2 vector variables u^^^ , . . . u^") , so the existence of critical points 
is given by the vanishing of a discriminant, which is a polynomial in the 
coefficients of the function and therefore in A. This polynomial function 
of A is the generalisation of the characteristic polynomial for a multilinear 
form. However, the vanishing of the discriminant does not guarantee that 
the critical points will be real, which we require since their coordinates are 
the real and imaginary parts of the original complex coordinates. Thus our 
characteristic equation is a necessary condition, but not a sufficient one, for 
the existence of singular values of the hypermatrix. In our study of three- 
qubit states we shall see examples of solutions of the characteristic equation 
which fail to be singular values. 
To summarise: 

Theorem 3.2. Let a : C^^ x • • • x C*^" —>■ C be a multilinear function of n 
vector variables: 

«(uW,...,u(")) = a.,...„«;^..K;:). 

For A e IR, define the real polynomial function a{X) : R^'^^ x • • • x R^'^" R 
(identifying with R^'^) by 

a(A) (u(3), . . . , u(")) = det[A^A - ||u(3)f . • . ||u(")f /] 
where the di x ^2 matrix A — A{u^^\ . . . , u^")) is defined by 

A - ■ — n- ■ ■ ■ 

Suppose A 7^ 0. Then the equations 

a,,..,j^..^...v!C^=M^ (3.8) 

have a solution with all u^''^ non-zero if and only if a{X) has a real critical 
point. If this is so, A satisfies the polynomial equation 

Aa(A) = 0. (3.9) 
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4 Three-qubit states 



We will now apply the general theory of the previous section to the simplest 
case after the familiar bipartite (matrix) case, namely the case of three qubits. 
Even here the characteristic equation of a general 2x2x2 hypermatrix is 
dauntingly complicated. Instead of the most general three-qubit state, we 
will consider a state in the generalised Schmidt form 



with a, b, c, d real and positive, a > b > c > d, which has been shown [2] to 
be a canonical form for three-qubit states. 

This should be trivial. The first step in finding the canonical form for 
any state |^) is to find its injective tensor norm gd"^)); this gives the coeffi- 
cient a. The remaining real coefficients b, c, d are found by a further process 
of finding extrema of the same function whose maximum is gf(|^')). Thus 
we might expect that if \^) is already in the form (4.1), its characteristic 
equation should simply have the solutions a, b, c, d; the problem should be 
like diagonalising a matrix which is already diagonal. Interestingly, it is not 
so simple; so much so that we will need to simplify further by taking / = 0. 

In this case the function a{X) of Theorem 3.2 is a function of a single 
vector variable z e C^, which we will regard as a function of the four variables 
{zi, Z2,zi,Z2)- It is a homogeneous quadratic in both z and z and is therefore 
of the form 



where tijki = tjiki and tijki = tijik- Suppose cex vanishes together with its 
derivatives at z = a, z = b. Let Dzax{z) be the discriminant of ax regarded 
as a function of z alone; this function of z is singular if z = b (having the 
singularity z = a); thus Dz&x is a function of z which has a zero at z = b. 
We will now show that the derivatives of -D^cka also vanish at z = b. 

Since ax is a quadratic function of z, its discriminant is the determinant 
of the 2x2 matrix Uij = tijkiZkZi: 



I*) = a|000) + 6|011) + c|101) + d\lW) + /|111) 



(4.1) 



Oix{z, z) — tijklZiZjZkZi 




Hence 



dzk 




(4.2) 



Now 



dax 




dzm 
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so for m = 1, 2 the vector with components {vm)n — tmnpqbpbq is orthog- 
onal to a and therefore is a multiple of the vector with components eijUj: 

tmnpqbpbq Cuitnyd-f. 

By the symmetry of tmnpq, 

It follows that Cm is a scalar multiple of e^r^'r, so 

tmnpqbpbq — (TExnr^ns^r^s 



for some scalar a. Now (4.2) gives 

d[Dzax],,. , da 

— ^= (b) = tijkiaiajbi = — = 0. 

Thus DzOixi^) has a singularity at z = b, and therefore its z-discriminant 
vanishes. This gives us the characteristic equation for A. 

Now we can find the equation for the injective tensor norm of the three- 
qubit state (4.1). Thus we will apply Theorem 3.2 to the function 

Q;(x, y, z) = aijkXiVjZk 

where i,j, k = 0,1 and the only non-zero values of aijk are 

Oooo — o, Ooii = b, Oioi = c, Olio = d, am — /■ 

This notation is consistent with that of Theorem 3.1, with z = u^^^ being a 
single two-component vector. The 2x2 matrix A{z) is 



A{z) 



azo bzi 
czi dzQ + fzi J ' 



and the function ax is 

aA(z,z) = det[A(z)U(z) - A'zW]. 

We write this as 

ax{z,z) = F{zo,zi)zo + G{zo,zi)zqZi + H{zo,zi)zl. 

Then 

D,ax{zo,z^) = G^-AFH 
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which is a homogeneous quartic function of zq and zi: 

D,ax - a(A)^ + p{X)zlzi + 'y{X)zlzl + S{X)zozl + e(A)^, 

in which the coefficients a, . . . ,e are all quadratic functions of A^. The dis- 
criminant of this is the same as that of the associated function of one variable, 

P{x) = ax* + f3x^ + 70;^ + 6x + e, 

namely 

A = 27B'^ - 

where A — ae 1 , 

4 12' 

^ _ a^e aS^ + /3^e 7^ 
~ ~6 16 ^ ^ ~ 216' 

Thus the characteristic equation A = has degree 12 in A^. 

The equation simplifies considerably if we put / = 0. This special case 
has been investigated by Tamaryan et al. [10], who have found some in- 
teresting solutions (in addition to A = a, 6, c, d) with an elegant geometric 
interpretation. It is of interest to see how these solutions emerge in our 
approach. 

If / = the coefficients /9 and 5 vanish, so that D^axi^Q^'zi) becomes a 
quadratic form in (zq,z^), with coefficients 

a = ^ahcd{}? -o?){}? -d^), 
7 = [ol - 4a2)A^ + 4a3A^ - 8CT4, 
e = 4a6c(i(A'-62)(A2-c2), 

where is the /cth elementary symmetric function of a^, 6^, (? and ^. The 
discriminant of the quartic becomes 

A(A) = -^(f-4ae)^ 

= -\^a%^c^d^{\^ - a^){\^ - }?){\^ - ^){\^ - d^)Q{\f. 
Here (5(A) is a quartic in A^ which factorises completely: 

Q(A) = \\^S''\^ - L2][4(52 _ ^^^^)^2 _ _ 2ahcd)\ 
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where L and S are the symmetric functions of a,b,c,d introduced in [10]: 

= {ah + cd)(ac + bd){ad + be), 
S'^ = {s — a){s — b){s — c)(s — d), 
s = l(a + b + c + d). 

The discriminant A(A) has a zero of order 4 at A = 0, but this is not 
covered by Theorem 3.2. In fact, if A = the equations (2.1) become the 
equations for a singularity of the function q;(x, y,z), and the condition for 
a solution is the vanishing of the hyperdeterminant of a = (ajjfc)- But the 
known formula for a 2 x 2 x 2 hyperdeterminant ([4] p. 448) gives 

hdet(a) = Aabcd, 

so A = is not a singular value. 

This leaves as candidates for the injective tensor norm of the state 

1^) = a|000) + 6|011) + c|101) + d\llQ), 

i.e. the stationary values of the distance of the state from the set of separable 
pure states, or the singular values of the hypermatrix aijk, 

X = a,b,c,d,D = ^ and ^' = ^ (4-3) 

where 

L'"^ = — 2abcd — {cd — ab){bd — ac){bc — ad), 
S'^ ^S^- abed 

— ^{a + b + c + d){a + b — c — d){a — b + c — d){—a + b + c — d). 

These are the solutions found by Tamaryan et al. [10]. It is not assumed that 
S^, <S"^ or L'^ are non- negative. Tamaryan et al. point out that if S*^ > 0, the 
positive real numbers a, b, c, d satisfy quadrilateral inequalities (each of them 
is less than or equal to the sum of the others), and S is the area of the cyclic 
quadrilateral of which they are the side-lengths, while D is the diameter of 
its circumcircle. They also give a nice geometrical interpretation of S' and 
D': if 5"^ > 0, there is a self-intersecting cyclic quadrilateral ABCD with 
side-lengths a, b, c, d in which the vertices lie round the circle in the order 
A, B, D,C; the diameter of the circle is D', and S' can be interpreted in 
terms of signed areas. 
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In order to find the injective tensor norm we need to determine which of 
the solutions (4.3) is the largest. First we note that 

-4a'S' ( -a' + b' + + d' + ^4 >0 



— \ —a -\- -\- c -\-a i = — 

4 \ a / 4 

where we follow the notation of Tamaryan et al. in defining (and Vb, Tc, 
similarly) . Prom this we get the geometrically obvious result 



4^2 



(no side of a cyclic quadrilateral can be greater than the diameter of the 
circumcircle) , which also applies if a is replaced by b, c or d. Thus if the 
singular value L/2S is real, it is at least as great as any of the values a, b, c, d. 
Similarly, 



— ( -a^ + b^ + + d^ = ^ 

4 \ a J A 



L'^-4a^S'^^^ ( -a' + l)' + c' + d'-=^) = ^ > 0, 
from which it follows that 



The singular value L' /2S' can be real if S'"^ < 0, since L'^ can be negative, 
but then the above inequality shows that it cannot be greater than any of 
a, b, c or d. 

Finally, to compare the sizes of the singular values D — L/2S and D' — 
L' /2S', we note that 

D'' -D' = ^AiL' - 2abcd)S' - L\S' - abed)] = ^(L^ - 2S') 
and we use the factorisation [10] 

L^-(a''+b^+c''+d^)S^ = l(a^+b''-c^-d''){a^-b^+c''-d'')(a^-b^-c^+d^). 

Since a > b > c > d, the first two factors are non-negative; hence if 5"^ > 0, 

D'>D if and only if + > 6^ + 

With this ordering of a, 6, c and d, the conditions for positivity of -S'^ and 
S'^ are: 

^2 > <^ a<b + c + d, (4.4) 
^'2 > <^ a + d<b + c. (4.5) 
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Hence the greatest real solution of the characteristic equation, Amax (|^)), is 



Ams«(|*)) = < 



'a if a > b + c + d, 

D if a<b + c + d but a + d > b + c; 

D if a + d<b + c and + d"^ < b^ + c^; 

D' if a + d<b + chut a^ + d'^>b^ + 0^. 



However, this is not necessarily the injectivc tensor norm of |^). For A to 
be a singular value of Gijk, it is not sufficient that it satisfies the characteristic 
equation (3.9); the associated singular vector has to satisfy reality conditions. 
If A satisfies the characteristic equation, we can find a singularity of the 
function Dzax{z) at some vector z. Fixing this value of z, and considering 
det[A(z)M(z) — A^z"^z/] as a function of z, we are guaranteed to be able 
to find a singularity at some vector z; but this may not be the complex 
conjugate of z. 

If X — L/2S the function of z becomes 

^ ^ , , G/bcdj / 9 9 9 9 9 9\ /a \ 

D^ax{z) = (r„r^ - ^Tanr^va ZqZ^ + r^r^ z^) (4.6) 
whose non-zero singularities are given by 

TaTdzl-nrcZl = 0. 

Substituting for zi gives 

det[>l(z)M(z) - A^z^z/] = Kzl {raVdZQ - 2y^rJ^^r^ ZqZi + nVczl) 

where is a constant factor. This is singular at z where 

^TaTd zo - ^/r^cZi = 0. 

Thus zq/zi and 'zq/zi are both equal to \fr^r^Jr^d-, which must be real if 
they are to be complex conjugates, so r},rclrard must be non- negative. With 
the ordering a > b > c > d, only can be negative, so the condition for D 
to be a singular value of aijk is Va > 0, i.e. 

,2 2 ,2 26cci ^2 2^1 bed , ^ 

+ + > or < - + . (4.7) 

a 2 a 

We note that this cannot be satisfied if a > b + c + d, when D is not real, so 
(4.7) is sufficient to determine whether the geometric measure is a or D. 

If A = L'/2S' the function Dzax{z) is given by (4.6) with ra,---,rd re- 
placed by r^, . . . , r^, and its non-zero singularities z satisfy 
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but now substituting for zi gives 



det[A(z)t^(z) - A^ztz/] = K'zl (r'y^zl + 2^ r'j'^j'^z^z^ + r'^^zi) 



so 



If these are to be complex conjugates of each other, they must both be 
imaginary; hence t'^Jt'^t'^ must be non-positive if D' is to be a singular 
value of aijk. 

But this cannot be true if D' is the greatest solution of the characteristic 
equation, for then a-\-d < 6 + c but a? -\-cP > b'^ + c^, so that ad <bc; hence 

r„ — —a +0 + c + d 



a 

/ 2 t2 I 2 I t2 2acd 
Tu = a — + c + d ; — 



/ 2 I j2 2 I 72 2a6(i 

r^ — a +b — c + d 

c 

2 I 1,2 I „2 j2 _ '^(J'bc 
d 



and r'a^ + b^ + c- d- 



<2d^- 


2bcd 
a 

2acd, 
b 


<o. 


(4.8) 


> 2c^ - 


>o. 


(4.9) 


>2b^- 


2abd 

c 


>o. 


(4.10) 


< 2o? - 


2abc 
d 


< 0. 


(4.11) 



Therefore D' cannot be the injective tensor norm of the state |^). 

Thus we reach the same conclusion as Tamaryan et al.: the injective 
tensor norm (?( I \1/)) of the normalised three-qubit state |^) = a|000)+6|011)-|- 
c|101) + (i|110) {a > b > c > d) is determined by the single quantity = 
1 - + 2bcd/a: 

If r„ < 0, then t/(|*)) = a; (4.12) 



^r.,o^^^.m)^o^,! l;^^;J^^^^;^ . (4.3) 



5 The singular vectors 

The singular vectors associated with the singular values a, b, c, d, D, D' (i.e. 
the solutions (u^^^, u*^^), u*^^)) = (u, v,w) of the equations (2.1), or equiva- 
lently the product states |u)|v)|w) having extremal overlap with |^)) are as 
follows. 

For the singular values a, b, c, d the results expected. In each case 

there is just one set of singular vectors (up to phase factors), corresponding 
to the three-qubit states |000), |011), |101), 110). For the other two singular 
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values D and D' the function Dzdix has two singularities, yielding either two 
sets of singular vectors or none. For D singular vectors exist if ra,rb,rc and 
are all non- negative; these are [10] 



u = 



V = 



W) 



4:SVab + cd 

y/r^\0) + \/r^|l) 
AS Vac + bd 

\/^|0) + V^g^dll) 

4:Sy/ad + be 



(5.1) 



gives (up to 



Changing the sign of any of the square roots 
phase factors) just one further set of vectors. 

For the singular value D' there are singular vectors if either one or three 
of r^, r'fj, r^, are negative, so that one of them has a different sign from the 
others. Suppose the odd one out is r^; then the singular vectors are given by 



|u) = 
|v> = 

|w) = 



y^|0) + v/^|i) 





ab 




-r' r' W) 

' a' c\ / 


AS^/bd - 


ac 




-r'ar'Jl) 


AS^Jbc - 


ad 



(5.2) 



Again, changing any of the square roots gives one further set of vectors. 



6 Conclusion 

We have shown, in principle, how to solve the nonlinear singular-value equa- 
tions which determine the generalised Schmidt coefficients of a pure multi- 
partite quantum state (the singular values of the hypermatrix of coefficients 
of the state in a product bases). The largest of these Schmidt coefficients 
is the injective tensor norm, closely related to the geometric measure of en- 
tanglement. We have shown that it satisfies a polynomial equation which 
generalises the singular-value equation of a matrix, and we have shown ex- 
plicitly how to find the characteristic polynomial for a three-qubit state (i.e. 
a 2 X 2 X 2 hypermatrix). We have solved the equation to find the singular 
values and singular vectors for a class of such states. 
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It is interesting to compare our solutions with the bipartite case. The 
singular values and singular column vectors of an m x n matrix A have the 
following properties: 

• The singular vectors and null vectors of A span C". 

• The singular values are all real. 

• Every solution of the characteristic equation is a singular value. 

• The singular vectors associated with a particular singular value A form 
a vector subspace of C", whose dimension is the multiphcity of A as a 
root of the characteristic polynomial. 

• Singular vectors associated with different singular values are orthogo- 
nal. 

The example we have studied shows that most of these properties do not 
extend to the multipartite case (i.e. to higher-dimensional hypermatrices) . 
The first property docs survive in our example; indeed, in some cases there 
are enough singular vectors that not only do the singular vectors |u), |v), |w) 
span the individual spaces, but their tensor products |u)|v)|w) span the 
tensor product space. This is related to the failure of the the orthogonality 
property of singular vectors associated with different singular values: not 
only are they not orthogonal, they may not even be independent. Thus 
there may be more singular vectors than in the bipartite case. In another 
sense, there may be fewer singular vectors: those associated with a particular 
singular value form a discrete set (apart from phase factors), rather than a 
vector subspace. This is a feature of the nonlinearity of the problem. 

Our example suggests that something may survive of the link between the 
multiplicity of a singular value (as a root of the characteristic polynomial) 
and the number of independent singular vectors. However, there can be roots 
of the characteristic polynomial which are not singular values; they may not 
be real, or they may be real but have no singular vectors associated with 
them. 
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